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We show that a single quantum emitter can efficiently couple to the tunable plasmons of a highly 
doped single-wall carbon nanotube (SWCNT). Plasmons in these quasi-one-dimensional carbon 
structures exhibit deep subwavelength confinement that pushes the coupling efficiency close to 100% 
over a very broad spectral range. This phenomenon takes place for distances and tube diameters 
comprising the nanometer and micrometer scales. In particular, we find a /3 factor 1 for QEs placed 
1 — 100 nm away from SWCNTs that are just a few nanometers in diameter, while the corresponding 
Purcell factor exceeds 10®. Our finding not only holds great potential for waveguide QED, in which 
an efficient interaction between emitters and cavity modes is pivotal, but it also provides a way of 
realizing quantum strong coupling between several emitters mediated by SWCNT plasmons, which 
can be controlled through the large electro-optical tunability of these excitations. 

PACS numbers: 78.67.Wj,73.20.Mf 


I. INTRODUCTION 

Achieving an efficient coupling between a single quan¬ 
tum emitter (QE) and the surface plasmons (SPs) sup¬ 
ported by metallic nanostructures has become a popular 
subject of research due to its potential application to 
quantum-optics [1-3] and sensing [4, 5]. This efficient 
coupling lies at the heart of several surface-based ultra¬ 
sensitive optical analysis techniques, which rely on the 
plasmon-driven enhancement of Raman scattering [4] and 
infrared absorption [5]. Remarkably, the localized SPs of 
a metal nanoparticle can enormously modify the spon¬ 
taneous decay of a neighboring excited molecule [6-8], 
while propagating SPs can produce similar effects over a 
broadband spectral range. Reducing the dimensionality 
of the plasmonic structure from 2D (metal surfaces) to 
ID (thin wires) enables better control over the coupling, 
which can be engineered to affect just a single SP [2]. 
Additionally, the SPs of ID geometries are well suited 
to act as mediators in the interaction between several 
QEs placed in close proximity to a plasmonic waveguide 
[9-12], thus suggesting the combination of these tools 
to design large-scale quantum-optics integrated devices, 
which could benefit from the plasmon robustness against 
environmental fluctuations to operate under ambient con¬ 
ditions. 

The recent emergence of graphene as a plasmonic mate¬ 
rial [13-15] has introduced an additional knob to improve 
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the performance of QE-SP coupling. The large electrical 
tunability and high degree of confinement recently mea¬ 
sured in graphene plasmons [16-20] has stimulated sugges¬ 
tions for their use in tunable plasmonic circuitry and meta¬ 
materials [21], as well as for the achievement of quantum 
strong coupling and efficient interaction with QEs [22, 23] 
with superior performance compared with conventional 
plasmonic metals. Although these plasmons have been so 
far observed only at mid-infrared and lower frequencies, 
their extension towards the more technologically appeal¬ 
ing spectral ranges of the visible and near-infrared has 
been argued to be attainable [24], particularly by reducing 
the size of the structures to scales of a few-nanometers, 
which are commensurable with existing graphene-related 
structures such as aromatic molecules [25] and carbon 
nanotubes. In particular, nanotubes of tens of nanometers 
in diameter have been recently suggested as suitable ele¬ 
ments for plasmon circuitry [26]. It should be noted that 
SWCNTs, like other carbon allotropes, exhibit UV plas¬ 
mons that have been well characterized in the past [27]. 
However, those plasmons are much lossier, and therefore 
less prone to efficiently couple to QEs, than the tun¬ 
able lower-energy plasmons on which we concentrate here, 
which only exist in doped structures and are predicted to 
display similar electrical tunability as graphene. 

In this paper, we show that quantum emitters can 
strongly couple to the electrically tunable plasmons of 
doped SWCNTs, reaching light-matter interaction levels 
that go even beyond those of planar graphene. The Purcell 
factor (i.e., the decay rate near the material, normalized 
to the decay rate in free space) is increased by nearly three 
orders of magnitude when reducing the dimensionality of 
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the carbon nanostructure from 2D (graphene) to quasi ID 
(nanotubes). More importantly, the coupling efficiency 
of the QE to the SPs supported by SWCNTs (i.e., the 
fraction of decay into plasmons, also known as (3 factor) 
reaches values nearing 100% over a very broad range of 
QE-SWCNT distances and QE/SP frequencies. 



FIG. SI. Sketch of the system under study. We consider 
a QE placed at a distance p — R from the surface of a SWCNT 
of radius R. Three orthogonal orientations of the emitter 
dipole p are considered, as shown by the red arrows. 


II. THEORETICAL FRAMEWORK 

Figure SI depicts the system under study: a single 
QE placed at a distance p from the axis of a SWCNT of 
radius R. The emission properties of the QE are deter¬ 
mined by its transition frequency uj and dipole moment 
p. We assume that R is sufficiently large as to neglect 
curvature effects (i.e., the discreteness of the electronic 
bands, as well as features near the Dirac point associ¬ 
ated with the tube chirality and finite radius). A recent 
study [24] indicates that this approximation works well 
to describe transversal plasmons for R > 1 nm under the 
doping conditions here considered (see below), although 
narrower nanotubes require a more fundamental level of 
description [28, 29]. SWCNTs have been synthesized in 
this size range [30] and their excitonic absorption bands 
well characterized [31-33]. We thus model the SWCNT 
as a hollow tube with the same surface conductivity aiuj) 
as graphene doped with the considered number of charge 
carriers per carbon atom. This allows us to readily com¬ 
pare 2D graphene and ID SWCNTs. For simplicity, we 
adopt the Drude model, cf{(jj) = {\e^E y // {oo + i/r), 
with realistic values of the Fermi energy Ep = 1 eV and 
the plasmon relaxation time r = 1 ps. Importantly, the 
relative comparison between both materials is indepen¬ 
dent of the choice of r. However, the high Purcell factor 
here predicted (see below) is roughly proportional to the 


assumed r. In this work we are most interested on situa¬ 
tions in which both R and distances p are small compared 
with the light wavelength A. Therefore, the electrostatic 
limit here outlined, which leads to relatively affordable 
analytical expressions, provides a very accurate level of 
description, as shown below. A comparison between re¬ 
sults obtained with the electrostatic approach and a full 
electromagnetic formalism [34] is presented in Fig. S2, 
whereas both approaches give nearly identical results on 
the scale of Figs. S3 and S4. 

We start our analysis by considering the screened in¬ 
teraction FF(r, r', cj), defined as the electric scalar po¬ 
tential created at the position r by an oscillating point 
charge exp(—icjt) placed at r'. Reciprocity implies that 
W (r, r', cj) is symmetric with respect to the exchange of r 
and r'. Also, it is convenient to decompose FF(r',r,cc;) = 
l/|r—r'|+FF^^^(r', r, uj) as the sum of bare and induced in¬ 
teractions. This quantity allows us to obtain the plasmon 
characteristics, as well as the decay rate of a neighboring 
QE. Direct solution of Poisson’s equation for both r and 
r' placed outside the tube yields (see Appendix) 


r', w) = - T (2 - <5mo) cos[m((y? - i^')] (SI) 

TT ^ 
m=0 

POO 

X / dk rm{k) cos[k{z - z')] Km{kp)Km{kp'), 

Jo 


where we use cylindrical coordinates r = (p, 2 ), (p), 

-ll{kR)Am 


T m{}J) — 


1 -h Im{kR)Km{kR)Ar^ 


(S2) 


is the reflection coefficient for cylindrical waves, = 
(47ri<j/ci;R)(m^ -h k^R^)^ and Im and Km are modified 
Bessel functions. The integral in Eq. (SI) is performed 
over the wave vector k parallel to the axis of the nanotube, 
while the sum runs over components of fixed azimuthal 
angular momentum number m. 


III. PLASMON DISPERSION RELATION 

Plasmon resonances are signaled by their strong re¬ 
sponse for a given external perturbation, or equivalently, 
by the poles of r^(/c). Here we should note that sign 
cancellations due to the exp (imp) modulation of the in¬ 
duced charge along the azimuthal direction of the tube 
surface render the contribution of m % 0 modes small if 
p is larger than the radius R. We thus concentrate on 
the dominant m = 0 plasmon band, whose complex wave 
vector /Cp is found as a function of frequency uj from the 
solution to the transcendental equation 

+ i/r) = Io{kpR)Ko{kpR) 

This dispersion relation agrees with previous studies that 
focus on the r ^ oo limit [35-37]. In Fig. S2(a) we show 
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Re{/Cp} as a function of light wavelength A for SWCNTs 
of radius R in the 2nni—100 /im range. This magnitude is 
normalized to the free-space light wave vector 
so that the plot directly illustrates the degree of spatial 
conhnement of the plasmon, whose radial electric held 
is proportional to K\{k^p\ Whereas the case of small 
radius corresponds to realistic SWCNTs, we also consider 
very large R in order to also deal with graphene coated 
cylinders. In this way, this hgure illustrates the evolution 
from small tubes to planar graphene with increasing R. 
For comparison, we also plot the dispersion relation of 
the SPs supported by a graphene sheet, clearly showing 
that plasmons propagating along carbon nanotubes are 
more conhned (larger Re{/Cp}) than those supported by 
graphene, and their conhnement increases with decreasing 
R. This suggests that SWCNTs are better suited to 
produce enhanced coupling with QEs. As discussed above, 
the electrostatic limit yields very accurate results for 
nanotubes of small radius {R < 100 nm), and also for 
larger radius at short wavelengths. 

The ratio between the real and imaginary parts of 
kp is also an important magnitude, typically used as a 
hgure of merit (FOM=Re{/i:p}/Im{A:p}) for evaluating the 
propagation characteristics of SPs. This FOM is plotted 
in Fig. S2(b) as a function of A for different SWCNTs of 
different radius R. Remarkably, these plasmons possess 
a larger FOM than those of graphene for spectral and 
geometrical-parameter ranges in which they also exhibit 
tighter conhnement, as noted above (see A < 100 jam 
region). These characteristics are very benehcial for the 
design of efficient coupling schemes between several QEs 
mediated by SPs. 


IV. PURCELL FACTOR 


A convenient way of assessing the strength of the QE- 
SP coupling consists in analyzing the rate of spontaneous 
emission, F. In particular, the Purcell factor P = F/Fq, 
where Fq is the decay rate in vacuum, is directly related 
to the ratio of the plasmon resonance quality factor to 
the mode volume. For a point dipole p located at r, it 
can be calculated as [38] 

where E^^^(r) is the held induced by the dipole at its own 
position, which can be in turn obtained from the screened 

interaction as E^^^(r) = — Vr[(p • Vr')lF^^^(r, r', cj)] 

r'=r 

Using Eq. (SI) and specifying for dipoles oriented along 
the three orthogonal directions shown in Eig. SI, we hnd 


Pp = ^ 


U m=0 


pOO 

W / dkk"^ [K'^{kp)f lm{rm{k)}, 

^_n J 0 

(S3a) 



FIG. S2. Propagation characteristics of plasmons sup¬ 
ported by SWCNTs and the limit towards graphene. 

(a) Real part of the plasmon wave vector kp as a function of 
the wavelength A of free-space light oscillating at the same 
frequency for different values of tube radius R. The plasmon 
wave vector is normalized to the light wave vector ko — 27r/A. 
The limit of 2D graphene (dashed curve) is smoothly approach 
at large Rs. (b) The corresponding hgure of merit (FOM) 
Re{/cp}/Im{/cp} of the guided plasmons studied in (a). Full 
electromagnetic theory (solid curves) is compared with the 
electrostatic limit (dotted curves) in both panels. 


) OO 

^0 Jo 


= 1 + / dk^Klikp) lm{rm{k)}, 


(S3b) 




P<fi = ^ + —Ts N Cm / dkKl^{kp) Im{rm{k)}, (S3c) 


where K^{z) = dKm{z)/dz, 6^ = (2 - 5mo), and = 

2w? I . 

In Fig. S3, we study the QE-SWCNT distance depen¬ 
dence of Pp, Pep, and Pz for tubes of radius R = 2 nm and 
for several QE emission wavelengths A in the 10 — 100 pm 
range. The highest Purcell factor is observed for radial 
orientation, as expected from the —1/p divergence of the 
induced held at small separations, whereas the azimuthal 
orientation renders poor coupling because it is only con¬ 
tributed by m 7^ 0 SPs. When comparing the Purcell 
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FIG. S3. Purcell factor. The main panel shows the Purcell 
factor [Eqs. (S3)] as a function of the distance p — R from 
the QE to the surface of a SWCNT of radius R = 2nm 
for several values of the free-space emission wavelength A 
and all three possible QE dipole orientations (see Eig. SI): 
radial (continuous curves), longitudinal (dashed curves), and 
azimuthal (dotted curves). The inset shows the corresponding 
Purcell factor in graphene for the same wavelengths and a 
dipole orientation perpendicular to the carbon plane. 

factor associated with either SWCNTs or graphene (see 
inset to Fig. S3) for the corresponding optimal QE-dipole 
orientations {Pp for SWCNTs and perpendicular to the 
2D carbon sheet for graphene), it is clear that not only 
the Purcell factor is higher for ID-SWCNTs than for 2D- 
graphene, but their spatial dependence is very different: 
whereas for graphene the Purcell factor increases rapidly 
towards short separations due to the dominant role of 
non-radiative channels, for SWCNTs this increase is much 
less pronounced. 


V. /3 FACTOR 

As mentioned above, a high coupling efficiency between 
QEs and propagating SPs is the key ingredient to achieve 
many of the proposed functionalities within the field of 
waveguide QED. In Figure 4 we show the (3 factor as a 
function of the QE-SWCNT distance for several values of 
the QE-frequency. This magnitude can be evaluated by 
calculating the contribution of the SP-pole to the integrals 
of Eqs. (S3). We evaluate the /3 factor near a SWCNT 
of radius R = 2 nm for the two optimal QE-orientations, 
radial (p) and longitudinal {z). For comparison, we plot 
in the inset the evolution of the coupling efficiency for a 
QE near a graphene layer, with its dipole oriented perpen¬ 
dicularly to the 2D carbon plane. The coupling efficiency 
for SWCNTs reaches 100% over a very broad range of dis¬ 
tances: for the optimal orientation (radial), P is close to 1 
for distances ranging from 1 nm to Ipm (spanning three 
orders of magnitude). Importantly, this high coupling 


efficiency extends over a very broad range of frequencies 
and is even larger in the very low frequency regime. Our 
results imply that quenching phenomenon (i.e., when the 
QE emission is dominated by non-radiative decay chan¬ 
nels) only appears at very short distances (< 1 nm) for 
QEs coupled to SWCNTs. As opposite to other metallic 
systems, SWCNTs present the crucial advantage of not 
requiring a spacer to avoid quenching of QEs at small 
separations. We attribute this negligible contribution of 
quenching modes to the very small quantity of material 
that binds the EM fields to the carbon nanotube. In 
contrast, in graphene a high coupling efficiency is only 
observed within a much narrower range of distances, and 
quenching shows up at distances below ^10 nm for the 
optimum wavelength, A = 10/im (see inset of Fig. 4). 



FIG. S4. [3 factor. Fraction of decay into plasmons {(3 factor) 
from radially (|| p) and longitudinally (|| 2 :) oriented QEs under 
the same conditions as in Fig. S3. The inset shows the [3 factor 
for a point dipole perpendicularly oriented to a 2D grapheme 
sheet. 


VI. CONCLUSION 

The strong interaction here predicted between QEs 
and the plasmons of SWCNTs opens new possibilities to 
implement waveguide QED schemes. It is worth noting 
that electrical contacts between gates and carbon nan¬ 
otubes have been extensively studied for their potential 
as nanoelectronics elements, thus facilitating the design 
of practical electrical doping schemes. This combination 
of classical electrical tunability and efficient coupling with 
QEs constitutes a powerful platform for the investigation 
of fundamental quantum physics and the design of devices 
capable of processing information encoded in the states 
of the QEs. 
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Appendix A: Quantum-mechanical versus classical calculations of the optical response 

Fig. SI demonstrates that our classical level of description works reasonably well down to relatively small radius 
and short light wavelengths. In this figure we compare the classical description for the nanotube, in which its dielectric 
response is modeled by a local conductivity a(uj) extracted from the local-RPA approximation (as in the main text), 
with a more microscopic theory. The latter is quantum-mechanical description that is based upon a tight-binding model 
for the electronic structure, used as input for the random-phase approximation [41]. In particular, the graph is obtained 
using a previously developed numerical code [42] for plasmons in ID graphitic structures with vanishing parallel wave 
vector. In both calculations the represented mode is the lowest-order one, which corresponds to transversal polarization 
of m = 1 azimuthal symmetry. In the main paper, we find that the dominant contribution to the decay of QEs to be 
produced by longitudinal modes (m = 0) integrated over a broad range of parallel wave vectors up to an effective cutoff 
^ 1/p, for which an even better level of accuracy is expected from the classical model, because the plasmon energies 
are then substantially below the Fermi energy, thus further away from the region in which coupling to electron-hole 
pairs can contribute to plasmon broadening and other nonlocal effects. 


Appendix B: Quasistatic approximation: Derivation of Eq. (1) of the main paper 

Electric and magnetic fields are decoupled in the quasi-static limit (c ^ oc), which is a reasonable approximation for 
structures such as the single-wall carbon nanotubes (SWCNTs) under consideration, where the dominant interaction 
between different regions of the material occurs at distances that are small compared with the light wavelength. In 
this limit, the magnetic field can be disregarded, while the electric field satisfies V x E = 0, so that it can be written 
as E = —V(/) in terms of an electrostatic potential (j). Using this expression. Coulomb’s law reduces to the Laplace 
equation = 0 in the empty regions inside and outside the SWCNT. The response of the material then enters 



FIG. S5. Comparison between quantum-mechanical and classical descriptions of the optical response of highly- 
doped SWCNTs. We represent the extinction cross-section normalized to the projected area of a (30, 30) SWCNT (radius 
2 nm) for light incidence and polarization directions both perpendicular to the tube axis. The Fermi energy is I eV. The 
quantum-mechanical description produces a broad, slightly blue-shifted feature compared with the classical model, which is a 
characteristic nonlocal effect and disappears for plasmon modes of energy well below the Fermi energy. 
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through the boundary conditions, namely, the continuity of the electric field components parallel to the tube, Ey, and 
the jump of 47rV • j/icj in the remaining perpendicular component right at the surface [43], where j = cr • Ey is the 
induced current and a is the 2D conductivity. 

We now present a succinct derivation of Eq. (1) of the main paper (the screened interaction from which 

the Purcell factor [Eq. (3)] readily follows from Eq. (2). We can calculate r', cj) as the potential induced at 

the position r by a point charge placed at r' and oscillating with frequency uj. We implicitly assume a e“^^^ time 
dependence in what follows and evaluate the conductivity a at that frequency. Given the symmetry of the structure, 
we use cylindrical coordinates r = (p, 2 ;, p), with the tube surface defined hy p = R. As we study the interaction of the 
SWCNT with dipoles located outside it, we only work out the solution for p, p' > R. 

We start by considering a complete basis set of solutions of the Laplace equation, Km{kp) and 

which are labeled by the real wave vector k along the 2 : direction and the azimuthal number m. The bare Coulomb 
interaction then appears as a projection over those solutions, 

^-- = -/ V (SI) 

where p< = min{p, p'} and p> = max{p, p'}. When is applied to the right-hand side of this equation, the jump at 
p = p' produces a 6{p—p') function, which can be readily combined with the closure relations Qirn{ip-cp ) _ 27r6{(p—ip') 
and f ^ = 27rS{z — z') to verify that this is indeed the Green function of the Laplace equation [i.e., 

V2(l/|r-r'|) = -47r6{T-T')]. 

In order to calculate the induced interaction we need to consider p> = p' for the point charge and p< = p near 
the tube surface, so we study the response to each component, which we regard as an external 

potential. Because of the axial symmetry, the corresponding induced potential must conserve k and m, and it must 
be Rm{kp) and outside and inside the tube, respectively. Any other combination of 

basis-set solutions produces divergences at either p = 0 or p —)► oc. These expressions implicitly define reflection and 
transmission coefficients Vm and which are determined by the boundary conditions at p = i?. In particular, the 
continuity of Ey also leads to a continuous potential, or equivalently, tmKm{kR) = while the jump in the 

radial electric field gives krmK'rni^^) ~ /uj){k^ ^ rn? jp^){tm + where the prime denotes 

differentiation with respect to the argument. The solution to these two linear equations produces 

-liikR)A^ 

™ l + Im{kR)Km{kR)Am' 

where Am = (dTria/u;+ k^R‘^). It should be noted that the derivation of this expression is facilitated by the 
use of the Wronskian Km{z)Im{z) — K'^{z)Im{z) z [44]. Einally, the induced interaction is obtained from the 
right-hand side of Eq. (SI) by replacing for /^, leading to 

y^md(r^ r', L 0 ) = - V] (2 - (5too) cos[m{tp - (f')] / dk r^ik) cos[k{z - 2 :')] Km{kp)Km{kp'), 

^ ^0 Jo 

[i.e., Eq. (1) of the main paper] after using the k ^ —k and m —m symmetries of the integrand. 


Appendix C: Results for a low doping level 

In this section we present similar results to those of the main text, but for a different doping level in the SWCNT. In 
Figures S2, S3 and S4, which are the counterparts for Figures 2, 3 and 4 of the main text, we use Ep = 0.5 eV instead, 
whereas the other parameters are exactly the same as those utilized in the figures of the main text. A conclusion of this 
analysis is that the capabilities of SWCNTs for ultra-efficient coupling to quantum emitters are very robust against 
changes in the doping level. Notice that in Fig. S4 we omit the calculation for A = 1.5pm, as there are no propagating 
surface plasmons supported by either graphene or SWCNTs for Ep = 0.5 eV at that wavelength (see Fig. S2). 
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FIG. S6. Propagation characteristics of plasmons supported by SWCNTs and the limit towards graphene, (a) 

Real part of the plasmon wave vector kp as a function of the wavelength A of free-space light oscillating at the same frequency 
for different values of tube radius R. The plasmon wave vector is normalized to the light wave vector ko = 27r/A. The limit of 
2D graphene (dashed curve) is smoothly approach at large R’s. (b) The corresponding figure of merit (FOM) Re{A:p}/Im{A:p} 
of the guided plasmons studied in (a). Full electromagnetic theory (solid curves) is compared with the electrostatic limit (dotted 
curves) in both panels. 


















Distance (nm) 


FIG. S7. Purcell factor. The main panel shows the Purcell factor as a function of the distance p — R from the QE to the 
surface of a SWCNT of radius R — 2 nm for several values of the free-space emission wavelength A and all three possible QE 
dipole orientations (see Eigure 1 in the main text): radial (continuous curves), longitudinal (dashed curves), and azimuthal 
(dotted curves). The inset shows the corresponding Purcell factor in graphene for the same wavelengths and a dipole orientation 
perpendicular to the carbon plane. 



EIG. S8. /3 factor. Traction of decay into plasmons (/3 factor) from radially (|| p) and longitudinally (|| z) oriented QEs under 
the same conditions as in Eig.SS. The inset shows the [3 factor for a point dipole perpendicularly oriented to a 2D graphene 
sheet. 
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Appendix D: Full electrodynamic calculation 
1. M and N functions 


The calculation of the electromagnetic properties in a system with cylindrical symmetry is facilitated by the 
introduction of two base functions: 

Mn,k,kM = V X[f{^u,], (SI) 

Nn,k,kA^ = — V X Mn,k,k,(f}, 

Qk 

where f{f) = Jn{kr) n ig an integer, and gk = . 

The N function corresponds to the TM polarization, while the M function is a TE mode. These functions satisfy 

V X Mn,k,k^ = 9kNn,k,k^ (S2) 

^ ^ ^n,k,kz ~ 9k ^n,k,kz 


and, in the coordinate system {ur^ue^Uz)^ they admit the following explicit expressions: 



(-Ukr)\ 

/ tkk, j:,{kr) \ 

^n,k,kz (^) — 

-kJ'n{kr) 

e^nee^k.z^ =- 

9k 

nk,, 

Jn\kr ) 
r 


\ 0 } 


\ k'^ J„{kr) J 


(S3) 


These wave functions are related to the even-odd functions found in the literature (for instance in Ref. [45]): 
M = Mg + iMo^ N = Ne -\- iNq. However the even-odd functions contain different angular dependences for different 
field components, and furthermore, they are plagued by changes of sign between even and odd components. The 
single-exponential functions are simpler, but contain the subtlety that simple conjugation does not provide an orthogonal 
set of functions. 

To perform the calculations we need orthogonal wave functions, which we find by inspection to be: 


^—inO ^—ikzZ 


^l,k,kS^ = [-y Jn{kr), -k J^{kr),0^ 

= ^ (-^kk,JUkr),-'^Ukr),eUkr)) e 

9k \ ' J 


(S4) 


—inO ^—ikzZ 


with the orthogonality relations: 


J ■ Nn,k,k. = 47r2fc 5{k - k') S{k, - k',) 5^ 

J drNl,k',k'S^ ■ Mn,k,kA^ = I rfr TVt, ,, ,, (rl ■ Mn,k,kM = 0 - 


(S5) 


At this point we find useful to simplify the notation: we define |T(o)), where o = {cr, n, /c, kz}^ and o = {<j, n, /c, kz}^ 
and a = such that 


(rivl/,) = Mn,k,kAr^. (4>o|r} = (S6) 

for cr = 1 = M, and a similar equation interchanging M and A for cr = 2 = A. 

The operators may act over different labels of T. For instance, the relations between M and A and the orthonormality 
expressions are written as 

= ^7^^k5{o - o') = NoS{o - o'), 

Vx |5'<,)=gfc|5'g), (S7) 

where Nq is an overlap factor (which could be included in the normalization of the functions M and A). 

We will also use the compact notation 


_ ^ pOO pOO 

do= E / / 


dkz • 


(S8) 


a=M,N n= — oo 
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2. Vacuum Green’s Function 


There are two dyadic Green’s functions commonly used, the magnetic and the electric ones, which are defined by: 

VX Vx^e-5^^e = l, (S9) 

V X V X Gm - Gm = V X 1, 
where g = uj/c\s the free space wavevector. 

It is more convenient to start by computing and derive Ge from it (Eq. 5.156 in [45]): 


Trm U'pU'p. 


Ge = ^ (v X G„ 

To compute Gm we represent it on the chosen basis set: 

G„ = y / dodo'ao,o'\'^o) (^'o 


(SIO) 


(Sll) 


Using Eq.(S8) and Eq.(SlO), we find — g‘^) Gm = V x 1, and projecting on (Ta,| from the left and |T^) from the 
right, we obtain: 

{^9kx ~ 9 ^ (^x,y — 9ky = 9ky^y^x,y 


Therefore 


9o 1 . 

^2 _ ^2 AT 


9o - 9^ No 


(S12) 


and 


Using Ge 


Gr, 




9o 


No - 9^ 


\^o) {^o\, 


(S13) 


-{urUr + V X Gm)l9^ and No = No, 9s = do, the expression for the electric dyadic Green’s function is: 


Ge = - 


UrpUrp 




9l 


5^ {gl - 5^) 


l^o) (^o|. 


(S14) 


An important point is that the integral is over waves with all values for k and even for those such that 
k‘^ k^ = gl g‘^^ which are not ’in-shell’. Now, the integral over k can be evaluated by contour integration in the 
complex /c-plane, which leads to: 


Jn{kr)Jn{kr') ^ i^f{k) f Jn{kr)H^^\kr'), if r' > r, 
k‘^ — {g‘^ — kl) 2k yH^^\kr')Jn{kr)^ if r'< r. 


(S15) 


where now k is ’in-shell’, i.e., it satisfies k^ k^ = g^. 


We now define an “in-shell” label 5 = {a, n, k^}, with f ds = X]cr=M n With this, we find 

= |: / ds^ 1^,) (^,-1, (S16) 

= 1 I f 1 

Gq = -^UryUry -|- J ds |4^s) (4^s| • 

Recall that these expressions come from an integral over /c, performed in the complex plane. When closing the 
contour, different contributions etc. appear depending on the relation between r and r'. Therefore, the 

expression (r|Ge|r') must be understood as: 


(r1Ge|U) = (rlT,)(T,|rO 


Ts(r) if r' > r, 

^'s (?05 if r'< r, 


(S17) 
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where is like ^s(r) by replacing J —)► 

We are interested on the dyadic Green’s functions for sources outside the nanotube, therefore the relevant vacuum 
Green’s functions are: 


GeO = -^Urtlr + -^ j ^ \ , (S18) 


3. Contribution to the Green’s function from reflected waves 


On top of the vacuum dyadic Green’s function, we must add a contribution in order to satisfy the boundary 
conditions at the SWCNT. In the region outside the SWCNT these contributions are associated to the reflexion of the 
"incident" wave |^s), which can be written as: 


Ge 


-i-JI 


L11 

with some coefficients rs^s' that must be determined form SWCNT properties. 
In the case of a SWCNT, rotational and translational symmetries imply 


(S19) 


s' — ^cr,cr' ^n,n' ^kz,k'^ 


and therefore 


Ger — - 

Stt 

KjTmr — ^ 


^ nOO -j 

E E / p I I ’ 

= -oo a,a'=M,N^-^ 

pOO -j 


(520) 

(521) 


gyp / ^ / V / 1,2 i^cr,n,l€z^ '^^cr',n,kz 

n= — oo a,a'=M,N " 

Inside the SWCNT, the fields should go as \^a,n,kz) (ftiey do not have to satisfy the radiation condition at r ^ oc). 
Taking into account rotational symmetry the expression for the dyadic Green’s function can be expressed in terms of 
the “transmission” amplitudes a' as: 


pOO -j 

y^ y^ / dkz-^tcr,a' \'^cr,n,kz) {^a']n,kz \ ^ 

= _nn rr rr'—A/T AT ~ 


Get = — 

Stt 

7 

_ ^ pOC -I 

^ EE/ \'^a,n,kj 


(S22) 


n= —oo a,a'=M,N 
oo 


n= —oo a,a' =M,N 

In order to find the coefficients rcr,a ', ta,a' we must impose the boundary conditions at r = i?. 


Ur X 
Ur X 


= 0 


(^1 - E 2 ) 

(#1 - # 2 ) = 2aE 


(S23) 


where a is link to the local conductivity of the SWCNT, a, a = 

The magnetic field is related to the electric field through H = —{i/g)V x E. On the other hand, from Eq.S19 and 
S22 we know that Gm = V x Ge- Thus, the magnetic field generated by a point electric dipole is: 


H = --Gm 
9 

and, thus the boundary condition in a SWCNT at the location of the carbon sheet is 

Ur X ^Gel — Ge2^ = 0 


(524) 

(525) 


X {^Gral Gfji2^ 


= 2ta Ge- 
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a. Reflection coefficients 


Incidence by a wave \M) 


— Ur X 


(GqI — Ge2^ 


= 0 


M 2 + = tNMN2 + tMMM2 

Ms + VmmM^^^ + = ^NmNs + tMIVtMs 

— Urp X {Gfril Gjyi2^!9 — ^Gjyi 

- A^3 + VmmN^^ + VnmM^^ - tNuMs - tMlVtNs = '^{tNMN2 + tMMM2) 
N 2 + + TnmM!^^ - t]SfMM2 - tMM^2 = ^(^ATM^S + 

where all components of N and M are to be evaluated at r = i?, and a = 2ia. 

In matrix form 


(S 26 ) 


(S 27 ) 




1 

to 

-N2 \ 

Mi') 

nP 

-M3 

-N3 

Tvi') 

Mi') 

—Ns T (y.M2 

—Ms T 0LN2 

Ui') 

Mi') 

— N2 — dMs 

— M2 — (xNs) 


/tmm\ 
^NM I _ 
tMM I 

X^NM / 


/ —M 2 

-Ms 

Ns 0 LM 2 
\—N2 + olMs/ 


• Incidence by a wave |A^) 

For this incidence a similar set of equations is obtained, with the same square matrix the vector of coefficients 
G = (^mat, ^ATAT, ^MAT, ^atat) , and a vector B with the same expression as before, but with the substitution 
N ^M. 

These systems of equations can be solved analytically. Defining adimensional quantities R = gR^ q = k/g and 
qz = kz/g^ we obtain: 


D = aHnJn{H'^Jn “ HnJn)^ Tz + 

qHKJn - HnJ'r, + aH„Jnq){aH'^J^ + q{H'^J^ - 

aJ^XlHKJn - HnJ'n + aHnJnq)R‘^ 


TmM = — - 


TnM = TmN = 


D 

(y.JnJn{ + HjiJ^nqqzR 

~D 

7 / ^3/^, TTf T' 


^NN 


aJl{{-H'^Jn + Hr,J'^)n^ql - q^aH'^J'^ + H'^J^q - H^R^)) 


D 


which, using the Wronskian relation: 




2i 


TTX 


simplifies to 


D = aHnJnWn^ql + q\W + aHnJnq){aH'^Jn + qW)!^ 
aX^q^iW + aHnJnq)R^ 


'f'MM — 


TNM = TmN 


D 

naRqqzW JnJ'n 


D 


tnn = - 


aJn\n^qlW + q^R‘^{aH^J^ + qW) 

D 


(S 28 ) 


(5 29 ) 

(5 30 ) 
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where W = Multiplying all numerators and denominator hy W we obtain: 


D = n^iraq^qRJnHn + q'^R^{l + iia q^RJnHn){l + na RJ'^H'^) 


(S31a) 


tmm = 


7raq^R^J!l^^{l ira q^RJ^Hn) 


D 


nira qzq^R^JnJ' 

2 (rr.2^2 , ^4 62 


'f'NM = ^MAT = — - 


tnn = — - 


iraqRJn" {n^ql + q^R^{l + TraRJ^H'J) 
D ’ 


where all Bessel functions have argument qR. 


4. Purcell factor 


In the CGS system, the spontaneous emission rate is 

r = ^ Im • Ge • p} (S32) 

A simple calculation shows that in vacuum 

r„„c = (S33) 

Therefore, the Purcell factor, defined as P = T /Tyac can be computed from the refiected part of the dyadic Green’s 
function as: 



E 





( 1 ) 


,n,kz 



(S34) 


As both rN,N and tm^m are even in n, but so are (^Tm) {^m\^ and (^Tat) The diagonal elements 

in N even in n. The non-diagonal elements vn.m and vm.n are odd in n, but the components of 

(TatI^ are also odd in n, so the partial contributions to the Purcell factor are all even in n. Similarly, using 
that, for all Bessel functions Jn{—x) = —Jn{x) and the fact that the Bessel functions always appear in pairs, the 
integral over kz can be restricted to the interval (0, oo). 

So finally, using adimensional units, we find 


o r oo 

P=l + -Re<^^(2-(5„,o) E 

j djQz 2 ip ' ^o',n,q^) i^a',n,q^ ' P) 1 

1 n=0 a,a'=M, 

nJo (1 J 


where 



d 

-qHyyqd) 


/ iqq^ (qd) \ 


, XN,n,qz = 

1 

^ 1 ^ 

1 


V 0 y 

^ q^Hi^\qd) J 


(S35) 


(S36) 


Recall that d = g with g = 27r/A, where d is the distance of the dipole to the SWCNT axis, and all dependence on 
the SWCNT radius appears in the refiection coefficients 
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5. Propagating surface plasmon modes supported by a SWCNT 


The dispersion relation of surface plasmon modes can be obtained from the poles of the reflection coefficients, given 
by the zeroes of the denominator D in Eqs.(S31). For small SWCNT radius, the relevant mode has n = 0. Therefore, 
N and M waves decouple and the pole appears in the coefficient, which is: 


r 


n=0 

NN 


ira q^R 


Jo{qR) 


2 


l + T:aq^RJo{qR)H^o^\qR) 


Therefore the condition for existence of SPP is: 

ql MqpR)H^^\qpR) = -C 

iraR 


(S37) 


(S38) 


Recall that qp = y 1 — For very confined modes, >> 1, so ^ tqpz is a complex number with a large 

positive imaginary part. Using Kn{z) = and In{z) = {7r/2)i^Jn{iz), we obtain the approximation, 

valid for c ^ oc. 


qlloiqpzR)Ko{qpzR) = (S39) 

ZQf R 

which coincides with the result obtained within the quasi-static approximation. 

In both cases, it is possible to write down the equations such that the left hand side depends on x = qpzR and the 
right hand side does not. For instance, using the quasi-static approximation. 


Io{x)Ko(x) 


icR 

A'KCT 


(S40) 


In order to estimate the plasmon wave vector, in the case of interest when the plasmon oscillates several times before 
it decays, we use the expression for a with r = oo (a = arriving at 


x'^ Io{x)Ko{x) 



(S41) 


where 



(S42) 


These expressions are simpler in both the limit for x = qpR « 1, when Io{x)Ko{x) ^ —log{x) and x = qpR » 1, 
when Iq{x)Kq{x) ^1/2. 
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